Exploring multipartite quantum correlations with the square of quantum discord 
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We explore the correlation distribution in multipartite quantum states based on the square of quantum dis- 
cord (SQD). In tripartite quantum systems, the necessary and sufficient condition for the SQD to satisfy the 
monogamy relation is given. Particularly, our analysis shows that the SQD can be monogamous for all three- 
qubit pure states, and we define a genuine tripartite quantum correlation measure which can be generalized to 
mixed states. Furthermore, we study the quantum correlation distributions in four-qubit pure states. Although 
these distributions are not always monogamous, they can be used to constitute correlation indicators. As a 
typical application, we use these indicators in the dynamical evolution of multipartite cavity-reservoir systems, 
which give an effective characterization of multipartite quantum correlations in the composite system. 
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I. INTRODUCTION 

Characterization of quantum correlation is an important but 
rather challenging issue in quantum information theory. As 
a typical kind of quantum correlation, quantum entanglement 
Jltl has been well understood in many aspects and widely ap- 
plied to quantum communication and quantum computation 
J2]. However, quantum entanglement can not account for all 
kinds of quantum correlations |3|]. It was shown that, using 
little or no entanglement, the deterministic quantum compu- 
tation with one qubit (DQC1) Jl can perform a task expo- 
nentially faster than any classical algorithm @|. There- 
fore, it is believed that the quantum speedup advantage comes 
from quantum correlation rather than quantum entanglement 
in the DQC1 model. Moreover, recent studies also identify 
that quantum correlation is a key resource in broadcasting of 
quantum states |7|,|8|], quantum state merging 1191, 1 1011 . assisted 
optimal state discrimination fiUl . remote quantum state prepa- 
ration II 1211 . and so on. 

Quantum discord (QD) is a prominent bipartite quantum 
correlation measure IU3l[l4ll . and has been accepted as a basic 
tool to characterize the quantum advantage beyond entangle- 
ment. Recently, generalization of the QD to multipartite sys- 
tems receives much attention. Using the relative entropy as 
a distance measure of correlation, Modi et al gave the defini- 
tion of QD for multipartite systems 11511 . Rulli and Sarandy 
proposed the global quantum discord by recasting the QD 
with relative entropies iflfSll . Moreover, with the difference 
of quantum versions of tripartite mutual information, two in- 
equivalent multipartite QDs were sugg ested id HI. How- 
ever, the characterization of quantum correlation structure in 
multipartite systems is still very challenging. Monogamy re- 
lation of entanglement is an important property in multipartite 
quantum systems, which means that entanglement can not be 
freely shared among many parties. As quantified by the square 
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of concurrences [19], entanglement is monogamous in multi- 
qubit systems 11201 12111 . and this property can be used to consti- 
tute genuine multipartite entanglement measures [20, 22, 23]. 
Therefore, it is natural to ask whether or not the quantum cor- 
relation is monogamous, especially for the QD. 

Streltsov et al proved that the monogamy relation does 
not hold in general for quantum correlation measures which 
are nonzero for separable states ll24ll . However, the quan- 
tum correlation measure can still be monogamous in special 
cases. For example, they showed that the geometric mea- 
sure of discord [25] is monogamous in three-qubit pure states 
ll24ll . The QD is a well-defined quantum correlation measure, 
and its monogamy property was investigated by many authors. 
Prabhu et al found that the QD does not satisfy the monogamy 
even for three-qubit W states [26]. Using the Koashi- Winter 
formula B27I1 . Giorgi pointed out that the condition for the QD 
being monogamous is equivalent to that for the entanglement 
of formation being monogamous l28ll . While Ren and Fan 
showed that, under the same measurement party, the QD is 
still not monogamous in three-qubit pure states 112911 . 

In this paper, we are motivated by the following two ques- 
tions: (i) whether or not the QD is monogamous in certain 
form? (ii) in what degree the discord is monogamous and can 
characterize the genuine multipartite quantum correlation? To 
answer the two questions, we explore the monogamy property 
of the square of quantum discord (SQD) in multipartite quan- 
tum systems. The paper is organized as follows. In Sec. II, we 
derive the necessary and sufficient condition for the SQD be- 
ing monogamous in tripartite quantum states. In three-qubit 
pure states, our results show that the SQD is monogamous, 
and we introduce a three-qubit correlation measure which can 
be further extended to mixed states. In Sec. Ill, we analyze 
the correlation distribution in four-qubit pure states and con- 
stitute a set of multipartite quantum correlation indicators. In 
Sec. IV, we use these indicators to characterize the dynamical 
evolution in the cavity-reservoir systems, which give an ef- 
fective characterization for multipartite quantum correlations. 
Finally, in Sec. V, we discuss some problems and give a brief 
conclusion. 



2 



II. MONOGAMY PROPERTY AND CORRELATION 
MEASURE IN TRIPARTITE QUANTUM STATES 

A. Definitions and monogamous condition 

In bipartite quantum systems pab, the total correlation 
can be quantified by quantum mutual information Ia-.b — 
S(pa) + S(p B ) - S(pab) with S(p) = -Trplogp being 
von Neumann entropy Jl3ll . While the classical correlation 
is given by J A -.b = max {E B } [S(p A ) - J2j Pj s (PA\Ef)l 
in which {E?} is a positive operator valued measure 
(POVM) performed on subsystem B and Pa\e b = 

-Yr B {EfpABE^)/ Pj with ft = TiAsiEf PAB Ef*) Q. 
The QD is used to characterize bipartite quantum correlation, 
which is defined as the difference between Ia-.b and Ja-.b, 
and can be expressed as fl3ll 

D A \B = S(pB)-S(pAB)+ttim {E B } ^2 Pj S(pA\Ef), (1) 

3 

where the minimum runs over all the POVMs, and Da\b is 
referred to as the discord of system AB with the measurement 
on subsystem B. The QD can also be written in the form of 
quantum conditional entropy Ir30ri 

D A \b = S(A\B) - S(A\B), (2) 

where the non-negative quantity S(A\B) = 
min^s j J2j PjS(pA\E B ) i s the measurement-induced quan- 
tum conditional entropy and S(A\B) = S(pab) — S(pa) is 
the direct quantum generalization of conditional entropy. In 
the following derivation, we will use the expression in Eq. (2) 
for simplicity. 

In tripartite pure states \i/j)abc, the entanglement of for- 
mation [19] and measurement-induced quantum conditional 
entropies are related by the Koashi-Winter formula 02711 

S(i\k) = S(j\k) = E f (ij), (3) 

where S(i\k) and S(j\k) are the conditional entropies 
with measurements on subsystem k, and Ef(ij) = 
min J2 £ PeS(Pi) is the entanglement of formation in the sub- 
system pij with the minimum taking over all the pure state 
decompositions {p e , p\^\ and i,j, k £ {A, B, C}. Using the 
Koashi-Winter formula, the quantum discord has the form 

A|fe = E f {ij) - S(i\k), (4) 

where the discord is under the measurement on subsystem k 
with k e {A, B, C}. 

Monogamy relation is an important property in multipar- 
tite quantum systems. As quantified by the square of con- 
currences, Coffman et al showed that the monogamy relation 
^a\bc — Cab+Cac * s satisfied in three-qubit quantum states 
and the residual entanglement can characterize the genuine 
tripartite entanglement H20I1 . However, for the quantum corre- 
lation, previous studies indicated that the QD is not monoga- 
mous even in three-qubit pure states 11261 l28l 12911 . 



Now we explore the monogamy property of the SQD. In a 
tripartite pure states \iI>)abc< we have the relation Da\bc = 
S{A) — Ef(A\BC) in which S{A) is von Neumann entropy 
and Ef (A\ BC) is the entanglement of formation under the bi- 
partite partition A\BC IU3U14I1 . Combining this relation with 
Eq. (4), we can obtain the quantum correlation distribution 
quantified by the SQD 

Da\bc ~ D a \b ~ D\\c 
= M(Ej) + 2S(A\B)[Ef(AC) - D(A\C)], (5) 

where M(Ej) = Ej(A\BC) - Ej(AB) - Ej(AC), and 
we use the relations S(A\C) = -S(A\B) and D{A\C) = 
E f (AB) + S(A\B) Q. In Eq. (5), the correlation distribu- 
tion contains two terms, in which the first is an entanglement 
distribution relation quantified by the square of entanglement 
of formation and the second is a function of S(A\B), Et(AC) 
and D(A\C). Therefore, the necessary and sufficient condi- 
tion for the SQD being monogamous is 

M(E 2 f ) + 2S(A\B)[E f (AC) - D(A\C)] > 0. (6) 

When the condition is violated, the quantum correlation dis- 
tribution is polygamous, i.e., D 2 A , BC < D 2 A , B + D 2 A > C . 



B. Monogamy property in three-qubit pure states 

In the following, we will focus on the quantum correlation 
distribution in three-qubit pure states. We first show that the 
term M(E 2 ) is nonnegative in the monogamous condition. In 
two-qubit quantum states, the entanglement of formation has 
an analytical expression Ej(pij) — h[(l + (1 — C 2 j) 1 ^ 2 )/2] 
in which h(x) — —xlog 2 x — (1 — sc)log 2 (l — x) is the binary 
entropy and Cy = max{0, \/Ai — \f\-2 — VA3 ~ V^} is 
the concurrence with the decreasing nonnegative A^s being the 
eigenvalues of matrix p^ ((J y ^a y ) p*j (a y ^)a y ). As a function 
of the square of concurrence, the entanglement of formation 
obeys the following relations 

E2 f( C A\Bc) ^ E f( C AB + C Ac) 

> E 2 (C 2 AB ) + E 2 (C 2 AC ), (7) 

where the monogamous relation C 2 A ^ BC > C AB + C 2 AC and 

the monotonic increase property of Ef(C 2 ) are used in the 
first equation, and we use the property that E 2 is a convexity 
function of C 2 in the second equation. According to Eq. (7), 
we can obtain M(E 2 ) > 0. 

Next, we analyze the second term in the monogamous con- 
dition via the generalized Schmidt decomposition. Under lo- 
cal unitary transformations, a generic three-qubit pure state 
can be written in the formyjj 

\^) abc = A |000) + A ie ^|100) + A 2 |101) + A 3 |110) 
+A 4 |1H), (8) 

where the real numbers Xi range in [0, 1] with the condition 
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FIG. 1: (Color online) The quantum correlation distribution 
of SQD in the generic three-qubit pure state as a func- 
tion of the parameters x(8i,8o) and 11(63, 82). The x- 
axis comes from parameters 80 and 81, and x(8i,8o) ranges 
in [x(0, 0), • ■ • , x(0, it/2), ■■■ , x(ir/40, tt/2), • ■ ■ , a?(ir/2, tt/2)]. 
The case for 1/(83, 82) is similar. The numerical results support that 
the SQD is monogamous in three-qubit pure states 



5^ A? = 1, and the relative phase <p changes in [0, tt]. Its two- 
qubit reduced density matrix of subsystem AB is 

/ 



PAB = 



A5 



A Aie~^ A A 3 \ 



A Aie^ Xl + X 2 2 g 

\ A A 3 g* A| + Xj ) 



(9) 



where the parameter g = XiXse 1 ^ + A2A4. The density matrix 
Pac is similar to pab, and we have the relation 



PAC = 5 



A 2 -H.A,- 



[P 



AB] 



(10) 



where 5a 2 «a 3 is a transformation interchanging the parame- 
ters A2 and A3. The analysis on the sign of the second term in 
the monogamous condition can be classified into three cases. 

Case 1. The parameters A2 = A3. In this case, we have 
the density matrices pab — Pac an d Pb = Pc- Therefore, 
the conditional entropy S(A\B) = S(AB)-S(B) = 5(C)- 
S(B) = 0, and the relation E f (AC)-D(A\C) = Ef(AC) - 
Ef(AB) — S(A\B) = 0. According to above relations, it is 
obvious that the second term in the monogamous equation is 
zero, and then the monogamous condition in Eq. (6) holds. 

Case 2. The parameters A2 7^ A3, but the following relation 
is satisfied 



In this case, we can derive that the second term in the monog- 
amous condition is non-negative, and then the SQD is monog- 
amous (see the proof in the part A of Appendix). 

Case 3. The parameters A2 7^ A3 and the relation in Eq. 
(11) does not hold. In this case, the second term in the monog- 
amous condition may be a negative value. Nevertheless, due 
to the first term M(E 2 ) is positive, it may offset the negative 
value of the second term, such that the SQD is still monoga- 
mous. 

The previous analysis on the monogamy property leads to 
a conjecture that the SQD is monogamous in three-qubit pure 
states. 

Numerical proof. All the three-qubit pure states can be writ- 
ten in the generic form in Eq. (8) up to local unitary transfor- 
mations. Therefore, it is enough to investigate the monogamy 
in the generic pure state. Without loss of generality, we can 
set Ao = cos6>o, Ai = sin6>ocos6*i, A2 = sin^osin^icos^, 
A3 = sin^osinfj'isin^cos^, and A4 = sin^osin^isin^sin^ 
with O , 9i,9 2 , 03 € [0, 7r/2] and 4> G [0, 7r]. We calculate nu- 
merically the quantum correlation distribution D 2 (A\BC) — 
D 2 (A\B) — D 2 (A\C) as a function of the five parameters 
Sq,Oi,02, 03 and </>, where the interval of is 7r/40 and the 
case for phase <fi is n/20. There are more than 2, 100,000 
quantum states in the numerical calculation and all the results 
is monogamous. In Fig.l, we plot the quantum correlation 
distribution as a function of 0q,0i, 02 and 03 (our numerical 
results show that the relative phase does not affect the corre- 
lation distribution and we set <\> = 0). The x-axis comes from 
parameters 0q and 9\, and the y-axis comes from parameters 
02 and 03. As shown in the figure, the distribution is zero when 
the parameters 0i — 0, 7r/2, which correspond to the tensor 
product states. For other cases, the distribution is positive. 
The numerical results support that the SQD is monogamous 
in three-qubit pure states. 



C. A genuine three-qubit quantum correlation measure 

A quantum correlation measure should satisfy the following 
necessary criteria: (i) it should be a non-negative real number; 
(ii) it is invariant under local unitary operations ll24l I32I . For 
an n-partite quantum correlation measure, it should still sat- 
isfy the condition: (iii) the measure is zero if the n-partite 
quantum state is product in any bipartite cut 0311 . 

Based on our previous analysis on the quantum correlation 
distribution of the SQD, we define a tripartite quantum corre- 
lation measure as 



Qai^ABc)) = D 2 AlBC 



D 



A\B 



D 



(12) 



\E f (AC)-E f (AB)\>\S(A\B)\. 



(11) 



which characterizes the genuine three-qubit quantum correla- 
tion in pure states. The nonnegative property of Q3 is sup- 
ported by the monogamy analysis on the SQD. The tripartite 
correlation Q3 is invariant under local unitary operations due 
to the QD being unchanged under the transformation. 

For the third requirement, we first consider the bipartite 
product state \ipABc) = \<Pa) ® \v>bc)- The SQD D 2 A ^ BC = 
S 2 (pa) = is due to the three-qubit pure state is a prod- 
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FIG. 2: (Color online) The genuine three-qubit quantum correlation 
Qs in the generalized three-qubit W state. 

uct state under this partition. The SQD E> 2 A \b = is be- 
cause we have (I a ® E?)pab{Ia ® ^f^) = Pab with 
Ej being the projectors composed of the eigenvectors of 
Pb- The case for D A , C = is similar. So, the gen- 
uine tripartite quantum correlation Q 3 (\ipABc)) = 0. For 
the product state \^'abc) = \VAb) <8> l^c), we a l so have 
Qs(Wabc)) = °- since ^a|bc = ^a|b = S 2 { PA ) and 
■^AIC = 0- Similarly, we can derive QsOV'abc)) = w i m 
I^abc) = \fAc) ® I^b). Therefore, Q3 is zero when the 
three-qubit pure state is a product in any bipartite cut, and then 
the third requirement for the correlation measure is satisfied. 

As an application of the correlation measure Q 3 , we con- 
sider the generalized GHZ state and the generalized W state, 
which are two inequivalent classes under stochastic local op- 
erations and classical communication 113411 . The generalized 
GHZ state can be written as \G 3 ) = a|000) + /3| 111). Its 
genuine tripartite quantum correlation is 

Qz{\G 3 )abc) = S 2 (A) = ha 2 log 2 a 2 -/3 2 log 2 /3 2 ] 2 , (13) 

where the correlation Q 3 is equivalent to the von Neumann 
entropy due to D A , B = D A , C = 0. For the generalized W 

state \W 3 )abc = o|001) + fojOlO) + c|100), its three-qubit 
quantum correlation is 

Q 3 (\W 3 )) = S 2 (A)-D AI b-D aic , (14) 

where S(A) = -c 2 log 2 c 2 - (1 - c 2 )log 2 (l - c 2 ) and = 
E f (ik) - S(i\j) mthi,j,ke {A,B,C}. 

In Fig. 2, we plot the correlation QsdW^)) as a function 
of parameters 8 and <f>, where we set the amplitudes a = 



cos6>sin0, b = cos(9cos0, and c = sin6> in the generalized 
W state. When the parameters 8 = 0, 7r/2 or <f> = 0, the gen- 
eralized W state is a product state and its genuine three-qubit 
quantum correlation Q 3 is zero. In other parameter regions, 
\W 3 ) is a nonproduct state and the correlation is positive. As 
shown in the figure, the maximal value is Q 3 {\W 3 )) ~ 0.2779 
which corresponds to the parameters 8 = <j> — tt/4. Our result 
shows that the W state has the genuine tripartite quantum cor- 
relation, although its corresponding tripartite entanglement is 
always zero. 



D. Correlation measure and indicator in mixed states 

In three-qubit mixed states, the quantum correlation distri- 
bution is not always nonnegative. As an example, we ana- 
lyze the quantum state Pabc(W) = \ipi){ipi\ + IV^XV^I, in 
which the non-normalized pure state components are \ipi) — 
a|100) + 6|010) + c|001) and \ip 2 ) = d|000), respectively. 
We set the parameters a = cos#i, b = sin#isin# 2 cos#3, 
b = sin6>isin0 2 sin6>3, and d = sin6'iCos6' 2 . When the param- 
eters 6\ = 62 = 8 3 = 0.47T, we can derive the correlation 
distribution D A , BC - D\, B - D\, c = -0.00517 (see the 
calculation in the part B of Appendix). 

Although the quantum correlation distribution can be neg- 
ative, we can still introduce a new measure for mixed three- 
qubit states via the convex roof extension of pure state corre- 
lation measure. The measure can be expressed as 

Qt(PABc) =^m{ Pi MPiQ3(lpABc)> ( 15 ) 

where the minimum runs over all pure state decompositions 
and Q 3 (iPabc) * s me measure in three-qubit pure states as 
shown in Eq. (12). According to the previous analysis on the 
pure state correlation measure, the mixed state measure can 
be nonnegative and quantify the genuine three-qubit quantum 
correlation. 

In general, the mixed state measure Q§ is difficult to calcu- 
late, since it involves not only the optimal pure state decompo- 
sition but also the optimal local measurement. However, for 
some specific quantum states, the measure is available. For 
example, we consider the mixed state 

PABC(P) =p\f)('p\ + (1 -jOI^XpH (16) 

where the components | ip) = (|000) + |lll))/v / 2and {p^-) = 
(|000) - |111))/V2 and the possibility p > 1/2. In a gen- 
eral pure state decomposition, its pure state component has the 
form \ip(q)) = y/q\f) ± ^/l~~q\p^) with q ranging in (0, 1). 
For this pure state component, we can get D AlBi^il)) — 
D A \cWq)) = and D A \BcWq)) = SA^(g)). Therefore, 
in the optimization of the mixed state correlation measure, we 
only need consider the optimal von Neumann entropy of qubit 
A. After some derivation, we can obtain that the mixed state 
con-elation is Q 3 (pabc(p)) = h 2 (l/2 + y/p(l - p)) with 
h(x) = — xlog 2 x — (1 — x)log 2 (l — x). 

Due to the difficulty in calculating Q§, we may introduce 
some alternative indicators for detecting the genuine tripartite 
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FIG. 3: (Color online) The quantum correlation indicators Q&(\ W4)) 
as a function of the parameters #1 and $2, in which the indicators 

(1*3 2) (2*2 2) 

Q\ ' and Q 4 ' do not work in the contour region of the bold 
blue circle but the indicator Q 4 2 * 2,3 ' works in all the parameter re- 
gion. 

quantum correlation in the mixed state pabc- These indica- 
tors have the forms 

Q 3 (A\BC) =max[0, Q 3 (A\BC)], 
Q 3 (B\AC) =max[0, Q 3 (B\AC)], 
Q 3 (C\AB) = max[0, Q 3 (C\AB)}, (17) 

where Q 3 (A\BC) = D 2 A ^ BC - D A , B - D 2 A ^ C is the quantum 
correlation distribution of mixed states in the partition A\BC, 
and Q 3 (B\AC) and Q 3 (C\AB) are the mixed state distri- 
butions in partitions B\AC and C\AB, respectively. When 
the indicator Q 3 is positive, it indicates that there is the gen- 
uine tripartite quantum correlation in the given mixed state. 
When the value is zero, the indicator does not work. It should 
be noted that when the indicator in partition A\BC does not 
work, the indicators in other partitions may work and detect 
the multiqubit correlation. For example, in the mixed state 
Pabc(W), we have shown that its quantum correlation dis- 
tribution in the partition A\BC is —0.00517 which results 
in that the indicator Q 3 (A\BC) does not work. However, 
when we consider the indicator in partition C\AB, we have 
Q 3 (C\AB) = 0.08406 which detects the multipartite correla- 
tion. 



III. MULTIPARTITE QUANTUM CORRELATION 
INDICATORS IN FOUR-QUBIT SYSTEMS 

In four-qubit pure states, the structure of quantum correla- 
tion distributions is more complicated than that in three-qubit 



states. We first consider the distribution 

Q 4 ^ = D A ^ BCD — D 2 A \ B — D A ^ C — D A \ D , (18) 

where the superscript number 1*3 means that the correla- 
tion distribution lies in the partition between one qubit and 
the other three qubits, the superscript number 2 means that 
the correlation of subsystems is in two-qubit forms, and the 
Dmbcd is referred to as the QD with the measurement on 
qubits BCD. 

(1*3 2) 

The correlation distribution Q\ ' ' is invariant under local 
unitary operations due to the quantum discords in the distribu- 
tion being unchanged under this kind of operations. More- 
over, for the product states \<Pa) ® \<Pbcd) and \ipab) ® 
Wcd), we can obtain Q^ 1 * 3 ' 2 -* = 0. In addition, under the 
qubit permutation, the correlation distribution is zero for all 
bipartite product states. For the generalized four-qubit GHZ 
state |G 4 ) = a|0000) + the con-elation distribu- 

tion is always nonnegative, and we have Q4* 3 ' 2 ' (IG4)) = 
(— a 2 log 2 a 2 — /3 2 log 2 /3 2 ) 2 . However, the case for a generic 
four-qubit pure state is different, and the correlation distri- 
bution may be polygamous which results in a negative value 
for the <2 4 1 * 3 ' 2 \ As an example, we consider the generalized 
four-qubit W state 

\Wi) = a|1000) + 6|0100) + c|0010) + d|0001> (19) 

where the amplitudes are set to a — cost^i, b = sin^icosf^ , 
c = sin^isini^cos-!^, andd = sin^isim^sim^, respectively. 
When #1 = 0.477T and $ 2 = #3 = 0.28-7T, we can derive that 
Q^' 3 ' 2 \\Wi)) = -2.1265 x 10" 4 (see the calculation in the 
part B of Appendix). 

According to different bipartite partitions in four-qubit pure 
states, we have other two quantum correlation distributions 
which can be written as 

^(2*2,2) _ j~)2 _ n 2 _ n 2 _ n 2 _ n 2 

«4 — U AB\CD u A\C 1J A\D ^ B\C B\D 

<2 4 ' - 1 = D 2 AB \ CD — D 2 AB \ C - D AB \ D , (20) 

where the superscript 2*2 means that the correlation distri- 
butions are considered in two-qubit partitions, and the super- 
script 3 means that subsystem correlations are in three-qubit 
states. Although the correlation distributions Q4S are not al- 
ways nonnegative, they can still be used to constitute multi- 
qubit quantum correlation indicators in four-qubit pure states. 
In this case, we define the following four-qubit quantum cor- 
relation indicators 

Qf^max^Qi 1 * 3 ' 2 )] 
Qi 2 * 2 - 2 >=max[0,Qf 2 ' 2) ] 
Q 4 2 * 2 < 3) =max[0,Q 4 2 * 2 ' 3) ]. (21) 

When the indicator Q4 is positive for a given quantum state, 
it means that there exists the genuine four-qubit quantum cor- 
relation. When the indicator is zero, it does not work and can 
not detect the multipartite correlation. 

For the generalized four-qubit W state in Eq. (19), we use 
the Q4S to detect the genuine four-qubit quantum correlations. 
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FIG. 4: (Color online) The indicators for multipartite quantum corre- 
lations in cavity-reservoir systems as a function of the time evolution 
nt and the initial state amplitude a. 



In Fig. 3, we plot the three indicators as a function of param- 
eters t?i and $2 (we set #3 = $2). In Fig. 3(a), the quan- 
tum states in the boundary correspond to the tensor product 

(2*2 3) 

states, and we find its four-qubit correlation Q\ ' ' is zero. 
The quantum states away from the boundary are non-product 
states, and we find that the indicator is always positive and 
detects the four-qubit correlation. In Fig. 3(b), the indicator 
Q\ is positive in most of the regions and can indicate the 
four-qubit correlation. In the contour region of the bold blue 
circle, the value of the indicator is zero, which means the indi- 
cator does not work for these quantum states. However, in the 

( 2*2 2) 

same region of Fig. 3(c), the indicator Q\ ' ' works. Sim- 

( 2*2 2) 

ilarly, when the indicator Q\ ' does not work (the con- 
tour region of bold blue circle), the indicators Q^* 3 ' 2 ^ and 
( 2*2 3) 

Q\ ' can detect the four-qubit correlation. In this sense, 
we argue that the four-qubit quantum correlation can not be 
characterized by one quantity and we needs a set of quantities 
to characterize it from different facets. Another example is the 
cluster state |C 4 ) = (|0000) + |0101) + |1010) - |HH))/2 
ll3~5ll . in which the four-qubit correlation indicator Q^ 1 * 3,2 '' = 
1 but the indicator Q 1 * 2 * 2 ' 2 ^ = 2. 

All the quantum correlation indicators Q4 are non-negative 
and invariant under local unitary operations. When the in- 
dicators are positive, they detect the four-qubit quantum cor- 
relation. Next, we will analyze the dynamical evolution of 
quantum correlation in four-partite cavity -reservoir systems, 
in which all the three-qubit and four-qubit indicators can char- 
acterize effectively the multipartite quantum correlations. 



IV. MULTIPARTITE QUANTUM CORRELATION 
CHARACTERIZATION IN CAVITY-RESERVOIR SYSTEMS 

Dynamical property of quantum correlations is very impor- 
tant in the practical quantum information processing, since 
quantum systems are always affected by the unwanted inter- 
action induced by the environment. Till now, the dynamical 
property of two-qubit quantum correlation has been widely 
investigated both theoretically and experimentally (see, for 
example, Refs. B36I439I1 and references therein). However, 
the dynamical property of multipartite quantum correlations 
is still very challenging. 

In this section, we consider the system composed of two 
entangled cavity photons being affected by the dissipation of 
two individual iV-mode reservoirs, where the interaction of a 
sing le cavity -reservoir system is described by the Hamiltonian 

N N 

H = huia^a + h~^^ujkb\bk + gk{ab\ + bk<^). (22) 

k=l k=l 

The initial state of the composite system is |$o) = (c*|00) + 
/3|ll)) ClC2 |00) rir2 , where the dissipative reservoirs are in the 
vacuum state. Under the limit of large N, the output state of 
the cavity-reservoir system has the form [40] 

|$t) - a|0000) ciriC2r2 + /0|&) ci n I0t)c 2 r 2 , (23) 

where \(f) t ) = £(i)|10) + x(t) |01) with the amplitudes being 
= exp(-«V 2 ) and X(t) = I 1 - exp(-Ki)] 1 / 2 . 
For the output state, we analyze its three- and four-partite 
quantum correlation indicators Q3 and Q4 given in Eqs. (17) 
and (21). In the calculation of other indicators, the multiqubit 
quantum discords can be obtained via the Koashi- Winter re- 
lation. However, it should be noted that the relation does not 
work for the calculation of some two-qubit quantum discord 
in the dynamical evolution, since it is involved in a solution 
to the entanglement of formation beyond the two-qubit case. 
Here, we use the method introduced by Chen et al for calcu- 
lating the quantum discord of X states (see the calculation in 
the part C of Appendix) 04111 . 

In Fig. 4, we plot the multipartite quantum correlation indi- 
cators as a function of the time evolution parameter nt and 
the initial state amplitude a. When the time Kt = 0, the 
quantum state is a product state and these indicators are zero. 
Along with the time evolution, these indicators are positive 
and can detect the genuine multiqubit quantum correlations. 
These indicators first increase to their maximum, and then de- 
cay asymptotically. When the parameter Kt — > 00, the output 
state evolves to a product state and all the multipartite quan- 
tum correlations disappear. 

In the cavity -reservoir system, its multipartite entanglement 
evolution was investigated in Refs. B40l 14211 . The genuine 
multipartite entanglement can be characterized by a series of 
entanglement indicators. Here, in our analysis, we consider 
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FIG. 5: (Color online) The indicators for multipartite quantum cor- 
relations as a function of the time evolution parameter Kt in com- 
parison to the multipartite entanglement indicators in the output state 
j$ t ) with the initial state parameter a — 1/vTO. 
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(24) 



where the subscripts i ^ j in the second equation. These mul- 
tipartite entanglement indicators are nonnegative and invariant 

(13 2) 

under local operations. The indicator E\' can be used to 
characterize the genuine multipartite entanglement in the par- 
tition c\\r\C2T2, and E^' 2 ' 2 ^ can indicate the genuine block- 
block entanglement in the partition c\T\ \c2r2 [42]. Moreover, 

(122) 

the indicator E$ ' is used to quantify the qubit-block en- 
tanglement in three-qubit mixed states IH3I445I1 . 

In Fig. 5, we plot the indicators of multipartite quantum cor- 
relations Q4 and Q3 as a function of the time evolution pa- 
rameter Kt in comparison to these multipartite entanglement 
indicators E4 and E3 in the output state |$ f ) with the ini- 
tial state amplitude a — 1/VT0. As shown in the figure, the 
multipartite quantum correlation is correlated with the mul- 
tipartite entanglement in every partition structure. However, 
the peak values of correlation and entanglement evolutions do 
not coincide completely. The reason is that quantum correla- 
tion and quantum entanglement are not equivalent in general. 
Particularly, in the dynamical procedure, the evolution of two- 
qubit entanglement can occur the phenomenon of entangle- 
ment sudden death [46-48], but the corresponding evolution 
of quantum correlation is always asymptotic. 



V. DISCUSSION AND CONCLUSION 



In a tripartite pure state \i/j)abc> the operational interpreta- 
tion of D A \ C is that the discord is equal to the total consump- 
tion in the extended state merging (ESM) from A to B i.e., 
the resource needed in the preparation the quantum state pab 
plus the consumption in the state merging 11 1 CHI . Our analy- 
sis and numerical results show that the quantum correlation 
distribution D 2 A ^ BC — D 2 A ^ B — D 2 A ^ C is monogamous in three- 
qubit pure states, which can be explained that the square of 
the entanglement in the partition A\BC is no less than the 
square sum of entanglement consumption in the ESM proce- 
dures from A to B and A to C. 

Unlike quantum entanglement, quantum correlation char- 
acterized by the quantum discord is asymmetric, which means 
that Da\ b ^ D b \a in general. In the quantum correlation dis- 
tribution D \\bc~ D 2 A \ B — D 2 A ^ C , we have analyzed the SQDs 
with the measurement on different parties. However, when we 
consider the measurement on the same party, the correlation 
distribution D 2 BC ^ A — B 2 B ^ A — D 2 C | A may give a different value 
and its property is worth to further investigate. 

In an A-qubit system, we can define different multipartite 
correlation indicators according to the structure of quantum 
correlation distribution. For the generalized A-qubit GHZ 
state \GHZ n ) = a\0)® n + we have an effective in- 

dicator 



i),2] _ n 2 



zZ D M\a^S 2 {A 1 ) 1 (25) 



which detects the A-qubit quantum correlation. For the gen- 
eralized A-qubit W state \W N ) = ai|10 . . . 0) + • • • + 
a k \0 . . . 1 . . . 0) H h a„|00 . . . 1), the indicator 



J(JV-2)*2,(JV-1)] 

'n 



= D 



A l ---A N _ 2 \A N - l A N 

An- 2 \Aiv 



D 



n 2 

U A 



A 1 ---A N - 2 \A N - 1 

(26) 



always works because the subsystem (A\ ■ ■ ■ Ajv-2) is equiv- 
alent to a logic qubit. 

In conclusion, we have explored multipartite quantum cor- 
relations with the SQD. In tripartite quantum systems, a nec- 
essary and sufficient condition for the SQD being monoga- 
mous has been given. Particularly, in the case of three-qubit 
quantum states, we have given the measure and indicator for 
the genuine tripartite quantum correlation. Furthermore, we 
have investigated the correlation distribution of SQD in four- 
qubit pure states and constituted a set of correlation indicators. 
As an application, these indicators are used to characterize the 
evolution of multipartite cavity-reservoir systems. Our work 
may be helpful for further understanding of the properties of 
quantum correlation in multipartite systems. 
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Appendix 
A. Proof of the condition in Eq. (11) 

We prove that the SQD is monogamous in three-qubit pure 
states if the condition in Eq. (11) is satisfied. The entan- 
glement of formation is monotonic increasing function of the 
square of the concurrence. When the entanglement in the sub- 
system AC is not less than that in the subsystem AB, we can 
derive the following relations 



E f( C Ac) > E i( C A.b) 



E f( C AB\c) t E f{C 2 AC\B) 

S{C) > S(B) 

S{A\B) > 0, (27) 



where we use the entanglement distributions C AB , C 



C 



AC 



C 



BC + t 3 and Ciq, B = C AB + C BC + r 3 with 
T3 being the three-tangle [ 20] . According to these relations, 
we know that the sign of Ef(AC) — Ef(AB) is the same as 
that of the conditional entropy S(A\B). Therefore, the second 
term in the monogamous condition can be written as 

2S(A\B)[E f (AC) - D(A\Q] 
= 2\S(A\B)\[\E f (AC) - Ef(AB) \ - \S(A\B)\].(2S) 

When the condition in Eq. (11) holds, this term is directly 
proportional to the absolute value of the conditional entropy 
S(A\B), which is obvious nonnegative. After considering the 
first term M(E 2 ) being also nonnegative, we get that the SQD 
obeys the monogamous relation 



D 2 (A\BC) - D 2 (A\B) - D 2 {A\C) > 0. 



(29) 



This completes the proof. 



B. Calculation of quantum correlation distributions 

In Sec. IID, via a three-qubit mixed state pabc{W), we 
show the quantum correlation distribution D 2 A ^ BC — D A ^ B — 

D 2 A \ C can be negative in the mixed state case. Now, we give 
the calculation. Based on the Koashi- Winter relation in Eq. 
(3), we have the discord D A{BC = E f {AE) - S{A\BC) in 
which the subsystem E is the environment degree of freedom 
which can purify the mixed state. Due to pabc{W) is a rank- 
2 quantum state, the environment subsystem is equivalent to a 
logic qubit. We can get the E f (AE) = 0.06942 by using the 
Wootters formula Then we have D\, BC = 0.10845. 

The discord for the qubits AB can be written in the form 
D A \ B = Ef{A\CE) - S(A\B), in which the entanglement 
Ef(A\CE) is available since the subsystem BE is equivalent 



to a logic qubit and then we obtain D 



2 

A\B 



0.02368. Sim- 



ilarly, we can derive D 2 A , C = 0.08994. Substituting these 
SQDs into the correlation distribution, it is found that the 
value of the distribution is —0.00517. 

In Sec. in, we show the quantum correlation distribution 
D a\bcd ~ D a\b ~ D A\c ~ D a\d can be negative in the 
generalized four-qubit W state. The calculation is as follows. 
The multiqubit discord D A \bcd is equivalent to the von Neu- 
mann entropy S(A) according to the property of pure state. 
By using the Koashi- Winter relation, we have the two qubit 
discord D A \ B = E f (A\CD) - S(A\B). Due to the subsys- 
tem CD being equivalent to a logic qubit, the Ef(A\CD) is 
available, and then we can deduce further D 2 A , B . The case 
for other two-qubit discords are similar. When the param- 
eters are chosen as #1 = 0.47-7T and #2 = $3 = 0.28-7T, 
the SQDs are D\ {BCD = 0.0053454, D A{B = 0.0018002, 
D 2 MC = 0.0018987, and D 2 A]D = 0.0018592, which results 



in the distribution Q y 4 



(1*3,2) 



-0.00021265. 



C. Calculation of the discord in cavity-reservoir systems 



The density matrix of two-qubit X state can be written as 



n AB - 

Px — 



( a 00 a Q3 
on ai2 
a* 12 a 2 2 



V 



'03 



a 33 



When the elements satisfy the following relations [41]: 

|ai2 + ao3| > K2 - ao3| 



IV a ooa33 - V a ll a 22| < |«12 1 



ao3 



(30) 



(31a) 
(31b) 



Chen et al proved that the optimal measurement for the quan- 
tum discord is a x . 

In the output state |$ t ) of the cavity -reservoir systems, the 
reduced density of two cavity photons has the form 



Pc l£ 



/ a 2 + f3 2 X * 

p 2 e x 2 



V a/3£ 2 



a^ 2 
,2,— > n I 






i ex z 
p 2 e 



It is obvious that the first relation is satisfied for the quantum 
state. Moreover, the second relation holds, due to 



K2I + I O03| - |V a 00«33 - V a H a 22| 

= V(go + gi) 2 - Jg$ + gl>0, (33) 

where go = a/3£ 2 and gi = fP^X 2 - Therefore, the optimal 
measurement for the quantum discord D C2 i Cl is er^ 1 , and the 
measurement operators are 



1 1 
1 1 



, Mi = - 
2 



1 -1 
-1 1 



(34) 



respectively. Then, according to the definition of the quantum 
discord in Eq. (2), we can get the value of E)^,, For other 
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two-qubit quantum discords in the quantum correlation indi- also a x . In a similar way, we can calculate these SQDs. 
cators Q 4 and Q 3 , we derive that the optimal measurement is 
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